stress.data <- read.delim("stress_data.csv")
attach(stress.data)
model.1 <- glm( psychological.stress ~ symptoms ); summary( model.1 )
model.2 <- glm( psychological.stress ~ income ); summary( model.2 )
cat("BIC for Model 1: ", BIC(model.1))
cat("BIC for Model 2: ", BIC(model.2))
cat("Support for model with smaller BIC relative to the other (K) =", abs(BIC(model.2)-BIC(model.1)))
